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I. INTRODUCTION 

The presence or absence of topologically protected 
states at the edge of a material is determined by the 
topology of its bulk Bloch states [T]. This topology 
can be characterized by integer invariants, that deter¬ 
mine whether the material is topologically trivial or non¬ 
trivial. Thus, development of methods to measure the 
topological invariants is one of the most important prob¬ 
lems in the field. The topological invariants of the elec¬ 
tron gas in external magnetic field (Chern numbers) are 
directly related to the quantized Hall conductivity [ 2 ]. 
Determining the topological invariants in various recently 
proposed counterparts of this system [3] is though more 
complicated. Particularly, the time-of-flight images were 
demonstrated to contain information about topological 
invariants of a cold-atom system ms]. Theoretically pre¬ 
dicted possibility to extract Zak phase, Chern numbers 
or Z 2 topological invariants from semiclassical dynam¬ 
ics of a wave packet HHia was recently experimentally 
realized for cold atoms in one-dimensional optical lat¬ 
tice m- Photonic systems are more preferable for re¬ 
alization of different measurement schemes due to eas¬ 
ier optical access to microscopic properties as compared 
to conventional electronic topological insulators or cold- 
atom systems. The methods to measure the topologi¬ 
cal numbers by tracing the individual edge states fin¬ 
gerprints in transmission spectra m or by manipulat¬ 
ing the single unit cell and directly measuring the Bloch 
function m were proposed for a lattice of coupled ring 
waveguides. The winding number of the scattering ma¬ 
trix eigenvalues was shown to determine the number of 
edge states and topological invariants pffilJ2Q] . although 
no clear way to measure them has been proposed yet. Re¬ 
cently, the relation between the surface impedance and 
Zak phase for a centrosymmetric one-dimensional pho¬ 
tonic crystal was revealed [21]. In this work we show 
that the topological indexes of a photonic system can be 
measured via the phase of the reflection coefficient, that 
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is accessible in simplest optical experiments. We consider 
both one-dimensional (ID) Aubry-Andre-Harper (AAH) 
photonic crystal [22H25] and two-dimensional (2D) lat¬ 
tice of tunneling-coupled resonators with synthetic mag¬ 
netic field [261ET], The direct correspondence between 
(i) the winding numbers of the reflection coefficient, (ii) 
the Chern numbers, and (iii) the presence of the edge 
states is proven. This provides the general recipe to ac¬ 
cess the edge states and topological indexes from outside 
the structure by optical means. 

II. TOPOLOGICAL INDEXES IN ID 
PHOTONIC CRYSTAL. 

First, we explain the proposed method for the ID pho¬ 
tonic crystal inspired by the AAH model [ 22 j [^. We 
consider the stack of alternating layers A and S, see 
Fig. la, where all the layers A have the same width 
while the width of the layers B is periodically modulated 
as [ 22 ] 

dB,n = <^23 [1 + 7 COS {27rbn + >c)] . (1) 

Here, is the width of unmodulated layers, 77 , 6 , and 
X are the modulation strength, frequency and phase; we 
focus on the rational frequency case b = P/Q^ where P 
and Q are the integers with no common factor. The 
difference between the dielectric constants of the lay¬ 
ers and sb leads to the Bragg reflection of the nor¬ 
mally propagating light. Namely, in the absence of mod¬ 
ulation the system forms an elementary photonic crys¬ 
tal with the period d = dj^ dB and the Bragg gaps 
in the energy spectrum around the integer multiples of 
ujb = {-s/ef\dAP-\/^dB)> The modulation (77 7 ^ 0) en¬ 

larges the unit cell size and the period from dto D = Qd, 
and drastically modifies the spectrum by splitting each 
Bragg gap into Q gaps. Importantly, the system exhibits 
cyclic evolution when the modulation phase x is continu¬ 
ously changed by 27r. In the real system, it might be pos¬ 
sible to realize this continuous variation if the modulation 
is induced by a running acoustic wave. This cyclic behav¬ 
ior allows one to map the ID system to 2D one, and intro¬ 
duce the topological indices — Chern numbers [23] [24] . 
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The considered ID photonic system is completely char¬ 
acterized by the dependence of the dielectric function 
on the coordinate along the growth axis z, and on 
the external parameter k. In the case of normal light in¬ 
cidence the wave equation at the frequency cj for the elec¬ 
tric field component E{z) perpendicular to 2; axis reads 

^E{z) + '^e^{z)E{z)=Q. ( 2 ) 

The properties of the eigenstates of Eq. 0 can be most 
conveniently analyzed in the reciprocal space. Perform¬ 
ing the spatial Fourier transform we obtain the system 
of coupled equations for the harmonics Ek-g with the 
wave vectors K — G, 

— K‘^)Ek + ^ cg^^Ek-g = 0 , ( 3 ) 

G 

where G = 2 'KglD (^f = 0, ±1,...) is the reciprocal lattice 

vector, q = u,/e^lc, Sg,^ = (1/-D) dze^(z)e~'^^, 

^G,>c = In the case ^ the quantity Cq ^ 

reduces to the structure factor, ^ oc exp(iG2:i), 
where Zi are the positions of the layers A inside the 
unit cell. The light with the wave vector q close to 
some reciprocal lattice vector G/2 exhibits Bragg diffrac¬ 
tion, leading to the formation of a stop band at the cor¬ 
responding frequency. In the spectral vicinity of this 
stop band we can use the two-wave approximation [ 28 ] 
valid in the regime of relatively weak spatial modula¬ 
tion of the dielectric function £^{z)^ which means either 
dA < dis^n or - eB\ < £aAb, so that |eG,><| < 1. 
We take into account in Eq. (§ only the two harmon¬ 
ics, E = {Ek^Ek-g)i with the wave vectors K ^ q and 
K — G ~ —q^ coupled by the structure factor component 
and neglect all other harmonics. This yields the 
equation 

dK,^ (T E=[q^-K^ + G{K - G/2)]E , (4) 

where 

dK,. = i-q^ Re -q^ Im G{K - G/2)], (5) 

and a is the vector of Pauli matrices. Introducing the 
dimensionless energy 5 = 4 g/G — 2^1 and the wave 
vector p = AKjG — 2^1, Eq. Q can be simplified to a 
Dirac-like form. 


P 


Ek 

= £ 

Ek 

-Eg -P _ 


Ek-g 


Ek-g 


where the structure factor plays the role of mass. 
Equation ^ yields the dispersion relation = p‘^ -\- 
describing the presence of the stop band of 
the width cG\eG centered at the frequency 

cG/(2V^). 

The topological properties of the problem are hidden 
in the dependencies on the modulation phase k. Partic¬ 
ularly, when K changes from 0 to 27 r, the system charac¬ 
teristics continuously vary and for k = 27 t they are the 


same as for x = 0 . Hence, the phase x can be treated 
as a wave vector in an auxiliary direction. This pro¬ 
vides the correspondence between the ID problem and 
the 2 D one [ 23 l [ 29 ], where the Bloch eigenstates de¬ 
pend on two wave vectors, K and x. Hence, we intro¬ 
duce the Chern numbers of the allowed bands in a stan- 
dard way as dx dK {dxA^ - d^Ax) /( 27 ri) , 

where Ak = / d2: ^£^0) £;*(K, x), 

E{K^ x) is the normalized solution of the wave equa¬ 
tion Q with the wave vector and is obtained 
replacing Dk by d^. Using the Dirac approximation §, 
valid in the vicinity of the stop band, we calculate this 
stop band contributions to the Chern number of the 
above- and the below-laying allowed bands. The results 
differ only by sign and are given by the solid angle swept 
out by dK^>c when K and x are varied divided by 47r [ 1 ]. 
This angle equals to the winding number of e^, i.e., the 
divided by 27 r phase that Cq gains when x is changed 
from 0 to 27 r. Thus, the allowed band is characterized 
by a Chern number given by the winding number of the 
structure factor in above-lying stop band minus that in 
the below-lying stop band. Conversely, the winding num¬ 
ber of the structure factor is a topological invariant of the 
stop band. Indeed, this number can not be changed with¬ 
out Cq ^ turning to zero for some x, which would mean 
eliminating the considered stop band. 

Now we show that the winding number of the structure 
factor Cq ^ coincides with that of the reflection coefficient 
Too from the semi-infinite structure for the incident light 
frequency lying inside the corresponding stop band. To 
this end the field inside the structure is expanded as a 
sum of Bloch waves and the field outside {z < 0 ) as a 
sum of incoming (oc e^^^) and reflected (oc waves. 

We apply the boundary conditions at the interface z = 0 
with the surrounding medium with the dielectric con¬ 
stant So and at z ^ oc. In the two-wave approximation 
Eq. ^ the resulting reflection coefficient is equal to the 
ratio of left- and right-going waves in the corresponding 
eigenstate of Eq. ?’oo(w) = ek-gIEk isniisi], where 
the wave vector K should cor respond to th e spatially de¬ 
caying eigenstate with p = - This relation 

directly links the measurable quantity, reflection coeffi¬ 
cient, to the topological properties of the Bloch states in 
the bulk. Using the Hamiltonian ^ we can obtain the 
reflection coefficient 


OO (^) — 


^G 




( 7 ) 


We see that for the energy £ lying inside the stop band 
we have full reflection, |roo| = 1, and the phase of the 
reflection coefficient is determined by the phase of struc¬ 
ture factor Cq Thus, the winding number of structure 
factor is equal to that of the reflection coefficient, and 
they both can be used to calculate the Chern numbers. 
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(b) Reflection coefficient intensity 



(c) Reflection coefficient phase 



FIG. 1. The color map of the squared amplitude (b) and the phase (c) of the reflection coefficient from the left edge of semi¬ 
infinite ID photonic crystal (a) as a function of light frequency u and the parameter x. Solid curves show the real part of 
the left-edge state frequency. The subplots in panel (c) show the cross sections of the phase map at different frequencies. The 
calculation parameters are b = 1/3, r/ = 0.5, n^/n^ = 2, dAlds = 0.2. 


III. PHASE SPECTROSCOPY OF 
TOPOLOGICAL STATES 

We have shown that from the phase of the reflection 
coefficient one can deduce the winding numbers of stop 
zones and, thus, the Chern numbers of allowed bands. 
Figure illustrates an application of this technique to 
the semi-infinite ID photonic crystal illuminated from 
the left. Panel (b) presents the absolute value of reflec¬ 
tivity, while panel (c) shows the phase of the amplitude 
reflection coefficient roo(<^)- Following the theory above, 
we characterize each stop band by a winding number i.e. 
the extra phase divided by 27r that the reflection coeffi¬ 
cient roo{uj) gains when the parameter x is varied from 
0 to 27r while uj remains inside the stop band. If we mea¬ 
sure the phase in some interval (e.g. from —it to tt, as 
in Fig. the winding number is equal to the number 
of phase cuts that cross the considered band gap, tak¬ 
ing into account their sign. Each cut connects two phase 
branching points. These phase branching points are the 
zeros of the reflection coefficient and can be located only 
in the allowed bands, since in the stop band |r(cc;)| = 1. 
Unless the gap closes, its winding number is invariant 
as there exists no way to eliminate a cut but annihilat¬ 
ing the corresponding branching points that is impossible 
while they are separated by the stop band. Comparing 
this result with the definition of the winding numbers 
through the phase of reflection spectra we find that the 
Chern number of a band equals to the number of phase 
branching points (zeros) of the reflection coefficient that 
are located inside this allowed band, taking into account 
their sign. 

Now we examine the bulk-boundary correspondence, 
i.e. we prove that the stop-band possesses edge states for 
any x, provided that the corresponding winding num¬ 
ber is nonzero. To this end, we analyze the analytical 


behavior of the winding number 

27r 

1 /'ain[roo(w + i7,xr)] 

«.(^ + .7) = 2- y - ^ -d- - (8) 

0 

in the complex plane of frequencies. Here, we have ex¬ 
pressed the winding number via the reflection coefficient 
from the semi-infinite structure and introduced the imag¬ 
inary part of the complex frequency 7. At the first stage 
of the proof we demonstrate that for every uj lying in¬ 
side the stop-band there exist some 7 > 0 and x such 
that roo{uj + i7, x) = 0 . Indeed, at 7 = 0 by defini¬ 
tion w{uj) 7^ 0 . On the other hand, for 7 ^ +oc the 
phase factors exp[(iu; — 'y)dA,B-s/£A,B/c], describing the 
propagation through the layers, are quenched and only 
the infinitely small front part of the structure contributes 
to the reflection coefficient. This part is ^-independent 
so the winding number vanishes. Since w{uj) 7^ 0, 
lim^^oo'^(<^ + 17) = 0, and the winding number is in¬ 
teger, we conclude that there exists a jump of the wind¬ 
ing number at certain 7. Such a discontinuity can only 
be caused by the singularity of the integrand in Eq. (§, 
that can be either a zero or a pole of the reflection coeffi¬ 
cient Too (<^ + 17, >c). The poles of the reflection coefficient 
for 7 > 0 are forbidden by the causality principle so the 
singularity corresponds to roo(<^ + 17, >z^) = 0. 

To finalize the proof we use the identity |roo(<^)P = 1 
valid for the real values of uj lying inside the stop band. 
Being analytically continued onto the entire complex 
plane it turns into roo{uJ + i7)r^(cc; — 17) = 1. Hence, 
the zero of the reflection coefficient ai uj-\-ij enforces the 
pole at UJ — 17. Such a pole means the presence of the 
radiative edge state with the frequency uj and decay rate 
7 [ 25 ]. The real part of the left-edge state frequency is 
shown in Eig. by solid curves. The points where the 
left-edge state appears or disappears as x is varied match 
the branching points of reflection coefficient phase. 
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All the above conclusions remain valid if we introduce 
some coating layers at the border of the photonic crystal 
with vacuum. Indeed, since the coating cannot close the 
stop band, the corresponding winding number remains 
unaffected and can be still used for calculation of Chern 
numbers and determination of presence of edge states. 
This reflects the topological protection of the considered 
states. 


IV. APPLICATION TO A 2D SYSTEM 


waves are given by the sum of the field over all the output 
waveguides at the corresponding edge, with the phases 
corresponding to the wave vector x, along the edge, which 
can be again realized by a proper tuning of the length of 
output waveguides. 

In the linear regime where {an,m) = Cin,m the system 
is described by the equation set 


+ 7)a 

n,m -v^ ( 9 ) 


Finally, we show how the phase spectroscopy can be 
used to measure the topological invariants in 2D systems. 
As described above, in the ID Aubry-Andre lattice the 
topological properties follow from the dependence of the 
reflection coefficient phase on the modulation parameter 
>c. Now we demonstrate that in a 2 D system the same 
topological invariants can be obtained by measuring the 
phase of the reflection coefficient from the edge of the 
system as a function of the wave vector along the edge. 

We consider a 2 D square optical lattice with one site 
per unit cell, see Fig. 2a. Its realization using the ring 
resonators linked by waveguides was studied theoretically 
and experimentally in Refs. [26| and |27| The nontrivial 
topology is induced by the synthetic magnetic field that 
stems from the specially engineered asymmetry of link 
waveguides that couple the neighboring sites. The direc¬ 
tion of the magnetic field is opposite for the clockwise 
and counter-clockwise resonator modes. For adiabatic 
links these modes are uncoupled, so in what follows we 
consider only one of them. The tight-binding Hamilto¬ 
nian for the given mode reads 

H — J ^ ^ l,m “1“ ) T H.C., 

n,m 

where J is the coupling constant, an,m and dj^ ^ are 
the photon annihilation and creation operators at the 
site and 27r6 is the phase acquired when hopping 

around the unit cell. We have chosen the gauge in such 
a way that for rational b = P/Q the magnetic unit cell 
has the shape Q x 1 . 

Figure shows the proposed scheme of measuring the 
topological numbers. We couple the probing waveguides 
to the sites on the left and right edges of a rectangular 
sample. The input waveguides on the left edge are all 
simultaneously coherently excited at the frequency uj in 
such a way that the phase difference between the adjacent 
sites is equal to x. As shown in the Fig. [^, this can 
be easily realized by using the same light source for all 
input waveguides and gradually changing the length of 
the input waveguide from site to site. Evidently, such 
an excitation scheme simulates an incident oblique plane 
wave with the wave vector x along the structure edge. 
The reflection and transmission signals are collected from 
the output waveguides at the left and the right edges of 
the structure, respectively. The reflected and transmitted 


where l^n,m — “1“ 

^n,N^) describes the light coupling to the probing waveg¬ 
uides, i.e. K,ri,m = if a probing waveguide is connected 
to the site (n, m) and Kn,m = 0 if not, is 

an incident field in the input waveguides on the left edge 
of the structure, that corresponds to the wave vector x 
along the vertical edge, and 7 takes into account the on¬ 
site losses. The reflection coefficient for the discussed 
geometry of detection reads 

Ny 

r(x) = ^ (1 + . ( 10 ) 

m=l 

The reflection coefficient intensity as a function of fre¬ 
quency and wave vector along the edge x is shown in 
Fig. One can distinguish two band gaps, where the 
reflectivity is close to unity. Fig. shows the phase of 
the reflection coefficient. Clearly, the winding numbers 
of the band gaps and the Chern numbers of the allowed 
bands follow from the phase map in Fig. in the exactly 
same way as for the ID system. 

Formal correspondence with the ID Aubry-Andre case 
is attained when the twist boundary conditions are intro¬ 
duced on the upper and lower edges of the system. Par¬ 
ticularly, we roll the rectangular lattice into a cylinder as 
shown by the dashed arrow in Fig. [^. The additional 
links with the twist phase — xA^, describing the flux of 
an effective magnetic field through the cylinder, are in¬ 
serted between the corresponding sites on the upper and 
lower edges. These twist links are described the extra 
terms Yln=i + H.c. in the Hamiltonian. 

Then the Fourier transform along the y axis shows that 
the only component excited is that with the wave vector 
X. This reduces the 2D system to a ID one, similar to 
one described in the first part of the work. Importantly, 
for large enough Ny the boundary conditions do not af¬ 
fect the result of the measurement. Thus we can omit the 
twist links and consider the unrolled rectangle geometry, 
as it was done for the calculation shown in Fig. ^,c. 

The multi-port scheme shown in Fig. can be trans¬ 
formed into a scheme with single input and output (see 
Fig. Il), which is expected to be easier in realization. 
The phase x that corresponds to the wave vector along 
the edge is now introduced by tuning the length of the in¬ 
put and output waveguides between the neighboring lat¬ 
tice sites in such a way that it corresponds to the phase 
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FIG. 2. Panels (a) and (d) show two proposed schemes of measuring of topological invariants. Panels (b) and (c) show the 
calculated color map of the squared amplitude and the phase of the reflection coefficient in the multi-port scheme for the 
rectangular x = 15 x 15 lattice, with b = 1/3, k,/J — ^jJ — 0.01. Panels (e) and (f) show the squared amplitude 

and the phase of the reflection coefficient in the single port scheme for the same parameters except for njJ — 0.01 and 
X iVy = 21 X 21. 


delay k for the input waveguide and —k for the out¬ 
put waveguide. The scheme is still described by Eq. 
where En,m = Sn,iEm + 5n,N^E!^ and 

Em = Em-ie'^ + (1 < TO < Ny) , (11) 

-E'm-1 =(l<TO<Vj^), (12) 
Eo = l, E'm = 0. (13) 

The reflection and transmission coefficients are equal to 
the fields and Eq, respectively. In the linear by 

K regime one can neglect the multiple tunneling of the 
light between the input/output waveguide and the lat¬ 
tice. Then the scheme presented in Fig. reduces 
to that shown in Fig. lit- The squared amplitude and 
the phase of the reflection coefficient in the single in¬ 
put/output configuration are shown in Figs[^ and f. One 
can see that the reflection phase in this configuration is 
similar to that shown in Fig. and can be also used to 
determine the topological indexes. 


V. SUMMARY 

To summarize, we have shown that the Chern number 
of a photonic structure can be deduced from the wind¬ 
ing number of the reflection phase. We have demon¬ 
strated that the nonzero winding number in a certain 
stop band gives rise to the topological edge states. In 
the ID Aubry-Andre-Harper photonic crystal the wind¬ 
ing number of the reflection coefficient in the stop band 
is equal to that of the structure factor characterizing the 
strength of the corresponding Bragg diffraction peak. To 
determine topological numbers of a 2D system one can 
use the reflection phase measured in the single- or multi- 
port configurations corresponding to oblique excitation. 
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